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We apply a quantum adiabatic evolution algorithm to a combinatorial optimization problem where 
the cost function depends entirely on the of the number of unit bits in a n-bit string (Hamming 
weight). The solution of the optimization problem is encoded as a ground state of the problem 
Hamiltonian H p for the ^-projection of a total spin-^. We show that tunneling barriers for the total 
spin can be completely suppressed during the algorithm if the initial Hamiltonian has its ground 
state extended in the space of the z-projections of the spin. This suppression takes place even if 
the cost function has deep and well separated local minima. We provide an intuitive picture for 
this effect and show that it guarantees the polynomial complexity of the algorithm in a very broad 
class of cost functions. We suggest a simple example of the Hamiltonian for the adiabatic evolution: 
H(j) = (1 — r) a* + t Hp, with parameter r slowly varying in time between and 1. We use WKB 
analysis for the large spin to estimate the minimum energy gap between the two lowest adiabatic 
eigenvalues of H(t). 

PACS numbers: 61.43.Fs,77.22.Ch,75.50.Lk 



I. INTRODUCTION 

Recently a novel paradigm was suggested for the design of quantum algorithms for solving combinatorial search and 
optimization problems based on quantum adiabatic evolution Q . In the quantum adiabatic evolution algorithm (QAA) 
a quantum state closely follows a ground state of a specially designed slowly varying in time control Hamiltonian. At 
the initial moment of time the control Hamiltonian has a simple form with the known ground state that is easy to 
prepare; at the final moment of time it coincides with the "problem" Hamiltonian Hp whose ground state encodes 
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the solution of the classical optimization problem in question. It can also be chosen to reflect the bit-structure and 
cost spectrum of the problem. For example, 



Here E z is a cost function defined on a set of 2™ binary strings z = {zi, . . . , z n } Zj = 0, 1, each containing n bits. The 
summation in (|l|) is over 2™ states |z) forming the computational basis of a quantum computer with n qubits. State 
\zj)j of the j-th qubit is an eigenstate of the Pauli matrix a z with eigenvalue 1 — 2zj = ±1). If at the end of the QAA 
the quantum state is sufficiently close to the ground state of Hp then the solution to the optimization problem can 
be retrieved by measurement. 

Running of the algorithm for several NP-complete problems has been simulated on a classical computer using a 
large number of randomly generated problem instances that are believed to be computationally hard for classical 
algorithms [|[ ||, || . Results of these numerical simulations for relatively small size of the problem instances ( n < 
20) suggest a quadratic scaling law of the run time of the quantum adiabatic algorithm with n. Furthermore, it was 
shown in j^] that the previous query complexity argument that lead to the exponential lower bound for unstructured 
search || cannot be used to rule out the polynomial time solution of NP-complete Satisfiability problem by QAA. 

On the other hand, a set of examples of the 3-Satisfiability problem has been recently constructed ||, [?], pi to test 
analytically the power of QAA in the situations where the optimization problem in question has multiple well-separated 
local minima. 

In these examples the cost function E z depends on a bit-string z with n bits, z = {z\, Z2, . . . , z n }, only via the 
Hamming weight of the string, w z — (z\ + Z2 + . . . + z n ), so that E z = f (w z ). The function f(w) is multi-modal, it has 
a local minimum separated from the global minimum by the barrier of an order- n width in w. For that reason classical 
local search like simulating annealing provably fails to find a globally optimal solution in time polynomial in n. It 
^ a "standard" QAA was applied to this problem in which a quantum evolution begins in a uniform superposition 
state — |z), ends in a target (solution) state |z t ), and the control Hamiltonian is a linear interpolation in time 
between the initial and final Hamiltonians. 

For the above examples it was shown || [jj that the system can be trapped during the QAA in the local minimum 
of a cost function for a time that grows exponentially in the problem size n. It was also shown [Q that an exponential 
delay time in QAA can be computed in terms of a quantum-mechanical tunneling for an auxiliary large spin system. 

It can also be inferred from fa, M that QAA will have an exponential complexity even if the cost function E z no 
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longer depends strictly on a Hamming weight but the deviation only occurs for states |z) that have exponentially 
small (in n) overlap with the adiabatic ground state wavefunction |^>o(i)) at all times during the algorithm execution. 

The above example has a significance more than just being a particular simplified case of the binary optimization 
problem with symmetized cost. Indeed one can argue that it shows one of the mechanisms for setting "locality traps" 
in the 3-Satisfiability problem But most importantly, this example demonstrates that exponential complexity of 
QAA results from a collective phenomenon in which transitions between the bit-configurations with low-lying energies 
can only occur by the simultaneous flipping of large clusters containing order-n bits. In many cases these transitions 
can be analyzed as a tunneling of spin variables. A similar phenomenon related to the tunneling of magnetization 
was recently observed in the large-spin molecular nanomagnets [[Io| . 

However low-energy collective behavior is also well known in spin glass models, many of which are in one-to- 
one correspondence with random NP-complete problems jnj. In particular, an important ingredient of the "replica 
symmetry breaking" picture of an infinite-range spin glass by Parizi Jl2| is that there are collective spin exitations 
that are of the order of the system size n whose energy is (0(1), i.e., it does not grow with the size of the system. A 



similar picture may be applicable to random Satisfiability problems |13|. 

Therefore in connection to the above example, it is important to understand how to design a polynomial time QAA 
without a prior knowledge of a particular form of the cost function / (w) (possibly multi-modal or even randomly 
sampled), so that a tunneling barrier between the local minima will be totally suppressed. This is a focus of the 
present paper 03. 



In Sec. Ill we present a theory of quantum adiabatic evolution of a large spin system introducing a control Hamil- 



tonian that guarantees polynomial time complexity for QAA in the symmetrized 3-Satisfiability example mentioned. 



In Sec. IV we provide an estimate of the minimum gap between the two lowest eigenvalues of the control Hamiltonian 



that determines the complexity of QAA. Sec. [v| contains concluding remarks. 

II. QUANTUM ADIABATIC EVOLUTION ALGORITHM 

In a standard QAA jjj one specifies the time-dependent Hamiltonian H(t) = H(t/T) 

H(t) = (1-t)H d +tH p , (2) 

where r — t/T is dimensionless "time". This Hamiltonian guides the quantum evolution of the state vector \ip(t)} 
according to the Schrodinger equation i d\ip(t))&t = H(t)\ip(t)) from t = to t = T, the run time of the algorithm (we 
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let h — 1). Hp is the "problem" Hamiltonian given in (|j). Hp, is a "driver" Hamiltonian, that is designed to cause 
the transitions between the eigenstates of Hp. In this algorithm one prepares the initial state of the system |^(0)) to 
be a ground state of H(0) = Hp. It is typically constructed assuming no knowledge of the solution of the classical 
optimization problem and related ground state of Hp. In the simplest case 

n-l 

H D = -CJ2H> W0))=2-"/ 2 ]T|z}, (3) 

j=0 z 

where <r J x is a Pauli matrix for j-th qubit and C > is some scaling constant. Consider instantaneous eigenstates 
|</)fc(r)) of H (r) with energies Ek(r) arranged in nondecreasing order at any value of r € (0, 1) 

H\<t>k)=Ek\<M, fc = 0,l,...,2 n -l. (4) 

Provided the value of T is large enough and there is a finite gap for all t £ (0,T) between the ground and exited 
state energies, AE(t) — E\{f) — Eq{t) > 0, quantum evolution is adiabatic and the state of the system \i/j(t)) stays 
close to an instantaneous ground state, \4>o(t/T)) (up to a phase factor). Because H(T) — Hp the final state \ij)(T)) 
is close to the ground state \4>q(t = 1)) of the problem Hamiltonian. Therefore a measurement performed on the 
quantum computer at t = T (r = 1) will find one of the solutions of combinatorial optimization problem with large 
probability. Quantum transition away from the adiabatic ground state occurs most likely in the vicinity of the point 
r w r c where the energy gap AE(t) reaches its minimum (avoided-crossing region). The probability of the transition, 
1 - \ {i)(t)\(j) Q {t/T))\ 2 t=T , is small provided that 

T> a"^2 » A£ min= mm [£i(t)-%)], (5) 

min U_ r _ 

(H T = dH/dr). The fraction in (|) gives an estimate for the required runtime of the algorithm and the task is to 
find its asymptotic behavior in the limit of large n> 1. The numerator in (|^) is less than the largest eigenvalue of 
H T = Hp — Hp,, typically polynomial in n However, A_E m j n can scale down exponentially with n and in such 
cases the runtime of quantum adiabatic algorithm will grow exponentially fast with the size of the input n. 

A. Binary optimization problems with symmetrized cost function 

Consider a cost function E z in the following form: 

n 

E z =.f(w z ), Wz^^Zj. (6) 
3=1 



This cost is symmetric with respect to the permutation of bits and w z is a Hamming weight of a string z. A particular 
example of this problem related to 3-Satisfiability was introduced in |5|, fj], || (the discussion in this subsection closely 
follows [0, Sec. 4). In this example 

E z = ^ c ( z * + z i + z <0> 

i<j<k 

c(m) = (l-8 m , )(l + 8 m , 1 (q-l)), m = 0,1,2,3 (7) 
where 6k,i is a Kronecker delta. For this particular case function f(w) in (^) takes the following form: 

f(w) — -^w(n — w)(n — w — 1) + —w(w — l)(n — w) + —w(w — l)(w — 2). (8) 
where q is an integer greater than or equal to 3. In the leading order inn> 1 one can write: 

/(«»)= Q)%(^)+0(n a ), (9) 

where 

h(u) = 4qu{l - uf + 4u 2 (l - u) + -u 3 . (10) 

3 

g(u) is a non-monotonic function with global minimum at u = corresponding to z±, z%, . . . , z n = Q. It also has a 
local minimum at u — 1 corresponding to z\, Z2, ■ ■ ■ , z n — 1 ( cf. Ref. Q), Fig. 1). 

In QAA the symmetrized cost function (||) corresponds to the following problem Hamiltonian Hp of total spin-n/2 
system (cf. Eq. (0) ) 

Hp = f(l- S z ) , (11) 
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where 



S 2 |z) =S z )|z), S, = --i«.. (12) 

Here S z is operator of ^-projection of a total spin § of a system of n individual spins 1/2. We used an obvious 
connection between the values of the Hamming weight function w z and corresponding eigenvalues S z of the operator 
S z . In what following we use "hat" notation for the operators of a total spin. 
It was shown H that if one uses Hu in a form corresponding to (|^) 

Hn=( n : 1 )a-S x ). (13) 



2 / V2 

then the minimum gap Ai5 min scales down exponentially with n implying the exponential complexity of QAA for this 
problem. 
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III. ADIABATIC EVOLUTION OF A LARGE SPIN 



A. Extended vs localized initial states 

In this paper we show that the main reason for the exponentially small minimum gap for the problem (Q) with Hp 



given in (13) is the fact that the ground state of the driving Hamiltonian is localized in the space of the z-projections 
of a total spin. We construct an example of the driving Hamiltonian with extended ground state and show that in 
this case the evolution time of QAA is polynomial in the number of qubits. 

In the following we adopt the notation for the total spin S = {^jli with j = x, y, z, where Sj are the projections 
of the total spin operator on the j-th axis. The total spin operator's components equal Sj = X]fc=i &j ano - are 
symmetric in all one-qubit spin operators {-S'j' '' j for k = 1, ...n. We use Sj for the eigenvalues of the operator Sj. 

We now introduce a new driver Hamiltonian 

H D = S 2 X . (14) 

Consider, for example, its ground state in the case when the total number of qubits n = 2s is even and therefore the 
total spin I = n/2 = s is an integer. Then the ground state is a state with S x = in a basis where x is chosen as 
a quantization axis. Making use of Wigner's rotation matrix d\ n m , (9) |20|, one easily projects this state onto the 
computational (problem) basis with the z quantization axis and m — S z . This gives us the ground state wave function 
in the problem basis as ijn) — d l m (tt/2). In the limit of large spin (which is the case of interest for us), the wave 
function is given by 



*o (m) « -juSi- m .2k, l=» is even). (15) 

VZn (l^ — m z ) ' L 

Note that ^ ( m ) = for the odd difference I — m = 2k + 1 and is given by (|5|) for / - m = 2k. From ([TJ), it follows 
that the ground state of S^ is delocalized in the m- space and spread over the whole interval, —l<m<l. 

On the other hand, the ground state of the operator (|l^ ) used in Q] as a driver, is localized in the m-space. Indeed, 
that ground state is a state with S x — I in the x basis. This gives us the ground state wave function in the problem 
basis, (m) = d l m ; (tt/2). In the limit of large spin one has 

1 / m 2 \ 

vEr' (m)«_ ex p(-- _), (16) 

which is clearly localized on the scale m 2 ~ /. The same conclusion holds for the case when n is odd and / take 
half-integer values. 



As we will show below, the adiabatic evolution of the delocalized states is fundamentally different from the localized 
ones. In particular, the localized ground states in general result in macroscopic tunneling. This gives the exponentially 
small in n ground state energy gaps and consequently the exponentially large complexity of QAA. Using the delocalized 
states one can avoid macroscopic tunneling. We argue that in a general situation when the information about the 
ground state is not used for constructing the driver Hamiltonian, the driver with extended ground state should lead to 
polynomial complexity of adiabatic algorithms independently of the specific form of the problem Hamiltonian provided 



that it is expressed as a function (in general, nonlinear) of the total spin operators, Sj (cf. Eq. (11)) 



B. WKB approximation for the large spin 



To be specific, we refer to the same problem Hamiltonian as in Eq. (|ll|). The full Hamiltonian takes the following 
form in the limit of large n 

H{T) = {l-T)nSl+T^)\(u), (17) 



n 



fi =ii-y, i=- 

2 I I ' 2 



where r = t/T ; < r < 1 and function h(u) is given in ([L0|). 

In order to get a simple physical picture, we will refer to the WKB-type approach commonly used in the theory 
of quantum spin tunneling in magnetics (QTM) jL7|, [jl8|, [lE|] . This approach is applicable for large spins {I » 1), 
which is the case of interest for us. We choose i as a quantization axis and following the standard procedure obtain 
the effective quasi-classical Hamiltonian in polar coordinates {9, ip} with 9 € [0, ir] and (p €E [0, 2ir]. In doing this, we 
make use of the following relations valid in the limit (I >> 1) 

5™ sa (^ sin sin </?)". 

As it was shown in [l9| , in the quasi-classical limit we have 

*— »-f. < 18 > 

which means that the motion of the quasi-classical spin can be described as a ID motion of a massive particle on a 
unit radius ring in the appropriate effective potential V (</>). We also have 

u=i(l-sin^). (19) 



Substituting this into (|l7|), we finally obtain the effective quasi-classical Hamiltonian in the form 

H(r) =n(l-r)p 2 +rg) 3 V(ip), 

V(<p) = h ( . (20) 



2 

One should note that the driver Ho — nS 2 does not change the effective potential caused by Hp and only introduces 
the effective kinetic energy into the problem. As we show below, this driver has extended eigenstates in the space of 
the z-projections to of a total spin, as opposed to the case Hd = n 2 S x considered in [jjj, where the eigenstates are 
localized in to. As we will see below, this leads to the absence of tunneling and polynomial gap in our case as opposed 
to the exponentially small tunneling amplitude arising in R . The adiabatic wave functions (f) satisfy 



^k (<p) = E k (r) fc (ip) , (21) 



where Ek (t) = Ek (r) /n is a rescaled energy. The first term in the l.h.s. of ( pl| ) corresponds to the driving 
Hamiltonian Hd and presents the kinetic energy of the particle moving in the effective potential due to the problem 
Hamiltonian Hp. At initial moment r = 0, the total Hamiltonian reduces to the driver Ho which describes a free 



massive particle on a ring and its eigenstates are well known ]20(]. In this case, the Schrodinger equation (|21[) gives 
exact wavefunctions and spectrum. Namely, we have 



*k(<p) = 4= ex P( l M, (22) 

V 27T 



Ek(0) = k\ 

where k is an integer number k — 0, ±1, ±2, ... for even n when the total spin I = n/2 is integer, and a half-integer 
number k = ±1/2, ±3/2, ... for an odd number of qubits n, when I = [n/2] + 1/2. One should note that in case of 
integer spin / , all eigenstates are twofold degenerate except for the ground state corresponding to k = and for 
the half-integer spin, all states including the ground state are twofold degenerate. This is a particular case of the 
Kramers' degeneracy pfl] , J18[ , which occurs due to the symmetry with respect to the spin-flip transformation. Note 



that since from (19) it follows that to = I simp, the states (|22j) are indeed delocalized in the TO-space. 

From (|2l]), it follows that the problem Hamiltonian is of the same order of magnitude as the driving one only for 
sufficiently short times rn 2 ~ 1. This means that the problem Hamiltonian is of the order of the level separation 
of the kinetic term only for sufficiently small times. Since the eigenstates of the driving Hamiltonian are extended 
(delocalized) in to, this implies the following qualitative picture of adiabatic evolution of the ground energy level with 
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FIG. 1: Multistable effective potential V vs polar angle tp for q = 3. 



r . At sufficiently small times rn 2 -C 1, the term due to Hp can be considered as a perturbation to the driver term 
and the energy levels are not strongly distorted. The eigenstates are delocalized in the m space. As r increases, the 
ground state is affected by the perturbation and the ground state energy increases. 



If the ground state is not degenerate (this is the case when the total number of qubits n = 2s is even), the gap 
may have a non-monotonic behavior in t in the range rn 2 ~ 1. Qualitatively, this can be described as follows. For 
sufficiently small times rn 2 -C 1, the ground state energy is increasing in r due to the diagonal matrix elements of the 
problem Hamiltonian until compensated by the level repulsion from the first excited state. After this, the ground state 
is pushed down and gradually approaches the ground state of Hp (which is E g = in the present case) . Since the 
upper levels are not as strongly affected by the perturbation as the ground state and since the number of upper levels 
is very large, O(n), the system of upper levels behaves as a rigid one. For this reason the strong interaction between 
the ground and first excited states occurs in the range of energies AE corresponding to the the level separation in 
the driver Hamiltonian and is linear in n. Therefore 



where AEp denotes the separation between the two lowest eigenvalues of Hp. Clearly, the energy scale of the minimal 
gap is polynomial in the number of qubits n provided that AEp is polynomial. 

As we have discussed above, in case of an odd number of qubits, all initial eigenstates of Hp, have a twofold 
degeneracy (Kramers' degeneracy). In this case, the degeneracy is removed by the effective potential (Zeeman splitting) 
and there is no avoided-crossing between the ground state and the first excited states of the total Hamiltonian H{f). 
One can easily show that in this case the level separation E\ (r) — Eq (t) grows monotonically in r [p0| (it is linear in 



C. Minimum gap analysis 




(23) 
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rj for rj <C 1). 

Apart from the features occurring in the range of small r ~ n~ 2 , the evolution of the ground state on the large time 
scale t ~ 1 is the same for both even and odd number of qubits. Most importantly, the estimate (^3|) holds globally 
in both cases. As we will see below, this is confirmed by numerical simulations. 



IV. MINIMAL GAP ESTIMATE 



As we discussed above, the ground state energy may have a non-monotonic time dependence on r when the 
total number of qubits is even and the ground state of Hjj is not degenerate. Because the eigenvalues of the driver 
Hamiltonian, n fc 2 , grow rapidly with the quantum number k it is possible to obtain a minimum gap estimate analyzing 
how Hp affects the two lowest eigenvalues of the driver. Making use of (0), we obtain the adiabatic gap as a function 
of time in the range rn 2 < 1 



(l-r)n-r(H P ) eu 



■8r 2 



(Hp) 



(24) 



Here we explicitely take into account that the first excited level of Hp, is twofold degenerate. Ssubscripts g, e\ t i above 
denote the ground state of Hp, and the two lowest exited states, respectively. Matrix elements of Hp on these states 
satisfy the following relations: 



( H p\, 9 = ( H r)e 



(Hp) g , ei = ( Hp )g,e 2 = ( H p) g , e 



(25) 



From (|24j), we obtain an estimate for the time when the minimal gap is achieved 



(Hp) 



9,e 



In the limit n»l we can write 



T) c = T c n 



( H P)e u e 2 



(26) 



(Hp) 



( H P)e u 



(27) 



Since matrix elements of Hp cc n 3 , it can be easily verified that the scaled quantity r\ c does not depend on n. The 
matrix elements can be calculated either in the quasiclassical basis given by (p2|) or exactly using (|l7j ) (see Appendix) . 
For the Hamiltonian (|l7j), the quasiclassical matrix elements yield 

9 « 



( H P)e u 



1024 



(28) 
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(for q = 3). Therefore, we have r? m i n = ^ ~ 7.1. This is in a qualitative agreement with ?7 m i n w 5.0 obtained in our 
numerical simulations. Substituting into we obtain the estimate for the minimal gap 

m 1/2 

AE min =AE(r c ) = nl-\ . (29) 

The corresponding value of the slope AE m - nl /n — s» 0.82 is again in a qualitative agreement with the value 0.86 
obtained in the numerical simulations (Fig. 5). 
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FIG. 2: Two lowest eigenvalues of H{t) @: i5o,i fs t 71^ in the vicinity of avoided crossing for n = 46 and g = 3. 
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FIG. 3: Difference between the two lowest eigenvalues of H(t) Ql_2) : E\ — Eo vs rn 2 n = 46 and q = 3. 



A. Numerical analysis 

We also performed numerical simulations of the adiabatic spectrum with Hamiltonian fllTf). In Fig. |[ we plot the 
ground state and the first excited state energies as functions of the dimensionless "time" r\ — rn 2 for even values of 
n. According to the above discussion, in this case the ground state is not degenerate and the energy level differnce 
exhibits a non-monotonic behavior (cf. Fig. |3|) leading to the polynomial minimal gap Ai? m i n ~ n at rj ~ 1. In Fig. ^ 
we plot an eigenvalue spectrum of the Hamiltonian H(t) at the avoided crossing point, r = r c for an even value of n. 
For k ^> 1 Ek — 0(n 3 k 2 ) which corresponds to the scaling analysis in Eqs. (|3l|) , (|3^) . Zeeman splitting of a doublet 
of the two lowest exited energy levels in Fig. H is of the order of n 3 r c = 0(n). 
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FIG. 4: Eigenvalues of the Hamiltonian H(t) (hj]) at the avoided-crossing, r = r c ; n = 46, g = 3. 



In Fig. 3 we plot the dependence of minimal gap AE m [ n vs n for even values of n . The insert to this plot corresponds 
to the dependence of the avoided crossing point r c vs n. 

Using the simple quasi-classical picture presented above it is not difficult to compute the minimum gap in the case 
of the driver Hamiltonian (13) in terms of the appropriate tunneling exponent and recover the answer given in 0. 
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FIG. 5: Minimum gap vs n for even n. The slope of the graph is 0.86 corresponding approximately to the analytical estimate for 
AEmi n = (2/3) 1//2 (p9|). Insert: Scaled position of the avoided crossing point r c n 2 vs n. The assymptotic result is r c n 2 w 4.98 
at n > 1. 
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V. CONCLUSION 

We show that macroscopic tunneling in QAA with the symmetrized cost function can be totally suppressed using 
a driver Hamiltonian with the ground state extended in the space of a total spin projection onto the direction 
of computational basis. This leads to a polynomial time complexity of QAA. We give a simple form of a driver 
Hamiltonian ( |l4| ) that has the aforementioned property and therefore makes the algorithm polynomial. 

We developed a simple intuitive picture of this phenomenon based on WKB approximation for large spins. It follows 
from the analysis that the results will hold for any form of f(w) that is sufficiently smooth on the scale w ~ 1. We 
conjecture that this phenomenon holds even if the cost function E z is symmetrical with respect to the permutation 
of bits but is not expressed only throught the Hamming weight of a string. 

We also argue using a picture of QAA as a quantum local search || |l6| that suppression of tunneling barriers with 
the operator changes the global neighborhood properties for QAA in a very profound way that can have an effect 
on the algorithm complexity for a larger class of cost functions, E z = f(w z ) + Af z where A/ z breaks the symmetry 
between the bits. In particular, it can be large for those states |z') that have exponentially small overlap |(z'|\I/o( T ))l 
with adiabatic ground states at all times (cf. Sec. 7.3 of ||). 

A possible generalization of the above analysis is related to random optimization problems with frustration, such 
as NP-hard problems and corresponding spin glass models. Exponential complexity of quantum adiabatic evolution 
algorithms for these problems is not necessarily related to tunneling but rather to the quantum diffusion phenomenon 
associated with the rapid falloff of correlations in the bit-structure with growing size of neighborhood around a given 
string p6[ | . The role of the tunneling and collective phenomena involving the low cost configurations in the performance 
of the quantum adiabatic evolution algorithms in random NP-hard is yet to be analyzed. 
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VII. 



NOTE ADDED 



Shortly after this work was completed, we learned about the work of E.Farhi, J. Goldstone, S. Gutmann |la| wherein 
another approach was suggested to achieve polynomial complexity of the quantum adiabatic evolution algorithm for 
the same optimization problem with the symmetrized cost function (ih. The algorithm proposed in Ref. ]l5fi is quite 
different from our algorithm. We believe that the proposal in [|l5| for random generation of interpolating "paths" H(t) 
in different trials of the algorithm is a very promising tool. It can perhaps be modified, by including intermediate 
measurents [^lj, to become an efficient adaptive algorithm for random NP-hard optimization problems. 

However we believe the particular method of achieving polynomial complexity presented in jl5| is less robust than 
ours in the problems with n qubits where the dynamics of the total spin n/2 is a key. The difference between the two 
algorithms is that in our approach a simple universal form of the driver Hamiltonian guarantees that the minimum gap 
scales polynomially in the problem size for a broad class of symmetrized cost functions f(w). Also, by construction 
our method does not require a specific knowledge of the solution, or even a specific form of the cost function f(w). 

We also believe that simple quasi-classical approach presented above enables one to study analytically a "volume" 
in the space of possible interpolating paths Jl5| ] that find a solution in polynomial time. Detailed study of this issue 
will be presented elsewhere. 

VIII. APPENDIX: MINIMUM GAP ESTIMATE USING MATRIX REPRESENTATION 

The eigenstates and eigenvectors of the Hamiltonian ( |l7| ) can be obtained solving from the matrix representation. 
Choosing a: as a quantization axis, we have S z = (—i) (S+ — S-) and 





(S + ) 



m,m — 1 



(S-), 



m — 1 , m 



-vW+I) 



m (m — 1)]. 



The full Hamiltonian (|l7| ) is given by 




and has the following matrix elements 
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(#) m , m -i = ir(|) 3 [(2-g)A 1 + ( 9 -2/3)A 2 ], (30) 

| / fi\ 3 
( H )rn, m -2 = 2^ (2 J (-9) A 3 , 

(^) m , ro -3 = ^(f) 3 (<Z-2/3)A 4 , 

with 



A (Z, m) = [Z (/ + !)- to 2 ] , 

Ai(J,m) = (/ + l)-m(m-l), 

A 2 a,m) = -^v^ (i + 1) - m ( m - 1) P (J + 1) - m ( m - 1) - 1] - 

Aa(I,m) = -^V[i (Z + 1) - m (m - 1)] [j (I + 1) - (to - 1) (m - 2)], 

A 4 (Z,m) = (Z + 1) - (to - 2) (to - 3)]A 3 (/, m) . 

One should note that the coefficients {A&} do not scale with I = 2n, i.e. Afe = 0(1) for n ^> 1. From (^0|), it follows 
that the eigenvalue problem H (r) ^>k = Ek (r) is rescaled for n ^> 11 in terms of dimensionless variables rj = rn 2 
and -Eft (r) = Ek (r) /n as 



H(r)*fc = £?fc(r)* fcl (31) 



with 



(£) = m 2 + ^(g+|-gA 

V / m,m lb \ 3 



16 V 3 

ft) = ^[(2-?)A 1 + (g-2/3)A a ] ) (32) 

m . m — 1 ID 



(h) = ^ (- q ) a 3 , 

V / m,m—2 lb 

(#) = ^(«-2/3)A4, 

\ / m.m— 3 ID 



implying that matrix elements of H do not scale with n and neither do the eigenvalues of H. It is strighforward to 
obtain an estimate for the minimum gap from the above equation using a Brillouin-Wigner perturbation theory in 
parameter to If one uses a 3-level truncation scheme, (to = 0, ±1), and sets q = 3 then a quasi-classical estimate (EE 
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can be recovered to the leading order in n. 
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